Theory and design of quantum cascade lasers in (111) n-type Si/SiGe 
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Although most work towards the reaUzation of group IV quantum cascade lasers (QCLs) has 
focused on valence band transitions, there are many desirable properties associated with the con- 
duction band. We show that the commonly cited shortcomings of n-type Si/SiGe heterostructures 
can be overcome by moving to the (111) growth direction. Specifically, a large band offset and low 
effective mass are achievable and subband degeneracy is preserved. We predict net gain up to lattice 
temperatures of 90 K in a bound-to-continuum QCL with a double-metal waveguide, and show that 
a Ge interdiffusion length of at least 8 A across interfaces is tolerable. 

PACS numbers: 73.43.Cd, 73.61. Cw, 78.45.-|-h, 78.67.Pt 
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I. INTRODUCTION 

Quantum cascade lasers (QCLs) have been developed 
in a variety of IH-V materials systems, although as yet 
there has been no successful demonstration in group 
IV. A Si/SiGc QCL would potentially reduce fabrica- 
tion costs and offer a route to photonic system-on-a-chip 
applications.'^ 

The most common approach towards a Si/SiGe 
QCL uses (001) oriented p-type structures^ and 
electroluminescence-'- '4 has been demonstrated. There 
are however significant challenges in designing p-type 
QCLs. The coexistence of heavy and light holes leads 
to fast nonradiative scattering, and strong valence band 
mixing causes large variations in transition energies and 
matrix elements with in-plane wave vector.^ 

The high longitudinal effective mass of SiGe A val- 
leys is commonly regarded as a major obstacle to n-type 
QCLs,^ and recent theoretical investigations have used 
transitions in the T and L valleys insteadi^i^ We have 
however shown previously that a (111) oriented Si/SiGe 
QCL using A valley transition is viable." 

In this paper, we compare the strain tensors for (001) 
and (111) oriented layers. We show that (111) oriented A 
valleys offer larger usable band offsets and lower quanti- 
zation effective mass than the (001) case and that compli- 
cations due to subband degeneracy splitting are avoided. 
We summarize our calculations of the principal scatter- 
ing mechanisms, current and gain, and present a bound- 
to-continuum QCL design in n-type (111) Si/SiGe. We 
calculate the waveguide losses and predict net gain in our 
design. Finally, we investigate the effects of temperature 
and nonabrupt interfaces. 



II. STRAIN TENSORS 

Lattice mismatch induces strain in thin Sii-^jGC;!; lay- 
ers on a relaxed Sii_a;^Gea;^ substrate, where x ^ Xs- 
Layers below their critical thickness deform elastically 
to match the in-plane lattice constant of the substrate 



and strain balancing of a multilayer structure is required 
to achieve mechanical stability. This is achieved by se- 
lecting a substrate alloy which minimizes elastic poten- 
tial energy with respect to in-plane straiufiS It is conve- 
nient to convert between both the interface coordinate 
system R — {x,y,z), where the z axis is normal to the 
layer interfaces, and the crystallographic coordinate sys- 
tem R' = {x' , y\ z'). For the (001) case, R' = R, whereas 
for (111) systems a transformation matrix U:R^ R' is 
required.JJ' 

The in-plane strain in R is defined as 



(1) 



where a is the lattice constant of the unstrained layer 
and fls is that of the substrate. A good approximation 
for lattice constant (in nm) ig^^ 



a{x) = 0.5431 + 0.01992a; + 0.0002733x2 



(2) 



Assuming isotropy over the xy plane, the strain tensors 
in R' for (001) and (111) oriented layers arei^ 
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where c-^ are the elastic constants, = — (c'^i + 2c'^2) 
and = c'li -I- 2c'i2 + 4C44. 

The minimum average strain energy corresponds to a 
substrate lattice constant 
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FIG. 1: Energy band schematic (not to scale) for a strained 
alloy layer on a rigid substrate. To find the potentials of the A 
valley minima in the strained layer, the difference in valence 
band energy and bandgap are calculated before adding the 
hydrostatic and uniaxial strain effects. 



where I is the layer thickness, k is the layer index and the 
elastic constants are grouped into a single term A. For 
the (001) and (111) orientations, the elastic constants are 
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III. BAND STRUCTURE 

As QCLs are intersubband devices, the indirect 
Sii-xGe^r bandgap is not an obstacle. The effective mass 
model for SiGe-based systems however, must account for 
transitions occurring away from the T symmetry point. 
The band structure is described as follows. 



A. Unstrained A valley potential 

The unstrained conduction band minima in Sii-a^Ge^; 
with X < 85% are located in six degenerate valleys in k 
space. Each valley lies close to an X symmetry point, 
and has a spheroidal equipotential surface with its major 
axis along the associated A directioufi^ 

The relaxed A valley potential relative to the substrate 
is determined from the model solid approximation.^^ The 
difference in average valence band maximum (in eV) is 



AVVb = (0.74 - 0.06a;^)(a; - x^),' 
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(8) 



and the highest valence band maximum is one third of 
the spin-orbit splitting above this, such that 

AVso 



AVvB = AVVb + 



(9) 
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FIG. 2: Usable range of A-valley offsets for (111) and (001) 
oriented Si wells with Sio.sGeo.s barriers as a function of Ge 
fraction in the substrate. For the (111) orientation the entire 
quantum well may be used whereas in (001) the Axy well 
potential defines the upper limit. 



where AVso is the difference in spin orbit splitting be- 
tween the two materials. 

The unstrained A valley conduction band offset is 
therefore 



AV( 
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where AEg is the difference in indirect F ^ A bandgap. 
To a good approximation, this is given (in eV) bj^ 



E„ = 1.155 - 0.43a; + 0.0206^2 



B. Strained layers 



(11) 



In general, two strain dependent terms are added to 
Eqn. [ini Firstly, hydrostatic strain, i. e. the overall vol- 
ume change, shifts the potential by an amount AVhyd = 
actr(e'), where Uc is the hydrostatic conduction band 
deformation potential. For (001) heterostructures, the 
hydrostatic shift is 
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whereas for (111) heterostructures, it is 
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Secondly, uniaxial strain, i.e. the distortion from the 
cubic lattice form, splits the in-plane A^j and the per- 
pendicular Az valley potentials by an amount 

AKni= (^±^)s„(43-£;i), (14) 
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TABLE I: Quantization and average two-dimensional density- 
of-states effective masses, and degeneracy n of A valleys in 
(001) and (111) SiGe alloys with Ge fraction under 85%, us- 
ing expressions derived from Ref. [Tsl . Masses are expressed 
relative to the rest mass of a free electron. 
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where 5u is the uniaxial deformation potentiaL The up- 
per sign represents the valleys and the lower sign rep- 
resents Axy This effect is absent in (111) heterostruc- 
tures as strain is identical along each of the principal 
crystallographic axes. 

The calculation is summarized in Fig. [1] The conduc- 
tion band edge is generally shifted by hydrostatic strain, 
and in the (001) orientation, uniaxial strain splits the A 
valley degeneracy. Although intervalley optical transi- 
tions are forbidden by the wave vector shift, degeneracy 
splitting complicates (001) QCL design by introducing 
nonradiative A^ A^y transitions. It is therefore desir- 
able to restrict optical transitions to energies below the 
Axy valley minima. 

For a strain balanced QCL with barriers narrower than 
wells, Xs < 25%. Fig. [2] shows that the usable energy 
range in the (001) orientation decreases almost linearly 
with Xs from a maximum of 150 meV at Xs=25%. The 
(111) orientation is therefore desirable for QCL designs 
as the entire A offset of 150 meV may be used regardless 
of the substrate composition. 



C. Effective mass 

The spheroidal A valley approximation remains valid 
in strained SiGe, as the effective mass only varies 
slightly. However, two separate effective masses are re- 
quired: the longitudinal effective mass rn;=0.916 for wave 
vectors along the major axis of a valley and the trans- 
verse effective mass mt=0.19 for wave vectors along the 
minor axes. The conduction band energy near a valley 
minimum is 

s2 / l2 fc2 ,2 \ 

EcB^\l^ + ^ + ^], (15) 
2 \mx my rriz J 

where rrix^y^z are constant effective masses for momentum 
in a given direction in i?, and k = k^x + kyij + k^z is the 
wave vector relative to the subband minimum. 

In QCLs, the quantization effective mass rUq describes 
the variation of conduction band potential in the growth 
direction, i.e. mq = m^. The density-of-states mass 
md accounts for in-plane motion and is anisotropic with 
respect to the in-plane wave vector. An isotropic approx- 



imation, md = ^m^my is commonly used howeverii2i2£ 
The effective mass values are summarized in Tbl. HI 

In the (GDI) orientation, the major axes of all four 
Axy valleys lie in-plane, so mq = mt and md = ^Jmtmi. 
Conversely, the major axis of each A^ valley points in 
the growth direction, so mq = m;, while is isotropic 
and equal to mt. 

It is slightly inaccurate to treat the A^ subbands as be- 
ing twofold degenerate in a simple effective mass approx- 
imation (EMA) as intervalley mixing splits the degener- 
acy in quantum confined systems. We have accounted 
for this in symmetric systems in a modified EMA and 
shown that the effect may be large in narrow quantum 
wells However, asymmetric structures require a com- 
putationally expensive atomistic calculation. ^^'^"^ As the 
splitting is at most a few meV in weakly confined states, 
the effect was neglected in the QCL simulations presented 
in this work. 

In (111) heterostructures, the situation is much sim- 
pler. All six A valleys have mq = 0.26 and identical 
confining potentials. Intervalley mixing is absent as each 
valley lies at a different value of k^y. The subbands are 
therefore sixfold degenerate within the EMA. 

In summary, the conduction band for (GDI) het- 
erostructures is complicated by effective mass anisotropy. 
The Az subbands may conceivably be used for QCL de- 
sign, but the large mq value severely limits the oscillator 
strength and hence the gain.— The A^^y subbands have 
lower mq, but states are weakly confined by the small 
band offsets and have low populations due to their high 
energy. 

The (111) orientation however, offers both low effective 
mass and high band offset and shows more promise for 
QCL designs. The low mass allows wider wells to be used 
than for A2 transitions in (GGl) systems. QCL designs 
in the (111) orientation are therefore more tolerant to 
deviations in layer thicknesses caused by growth errors. 



IV. CARRIER TRANSPORT MODEL 

Having established the differences in band structure 
and effective mass between (001) and (111) oriented het- 
erostructures, carrier transport may now be modeled 
in QCL structures. Detailed quantum theoretical ap- 
proaches have been used to simulate intersubband opti- 
cal emissions in quantum wells^^ and carrier transport in 
limited numbers of subbands in QCLsi^ They are, how- 
ever, too computationally demanding for use as design 
tools for large multi-level QCLs. 

Reasonably good agreement has been achieved be- 
tween experimental results and Boltzmann or rate equa- 
tion based models of bound-to-continuum THz QCLs in 
III-V systems i^'^1^^ In the present work, we have there- 
fore determined subband populations using a computa- 
tionally efficient rate equation approach. This has been 
described in more detail previously^ although a brief 
summary follows. 
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Several intravalley scattering mechanisms are impor- 
tant in Si/SiGe systems. Interface roughness scat- 
tering was calculated using the correlated Gaussian 
roughness model, modified for arbitrary interface 
geometries^^ and alloy disorder scattering was calcu- 
lated using a point perturbation modelj^*^ Ionized im- 
purity scattering was determined using a Coulombic in- 
teraction model as described by Unumai^^ with Thomas- 
Fermi screening)^ while electron-electron scattering was 
treated as a screened Coulombic interaction as described 
by Smeti^ The intravalley deformation potential scat- 
tering for electron-acoustic phonon interactions was also 
included. 

Intervalley electron-phonon scattering was determined 
only for the Si-Si branch of the deformation potential in- 
teraction, as the Ge fraction in quantum wells is small. 
The / processes, which transfer electrons to the perpen- 
dicular valleys, are faster than g processes which transfer 
electrons to the opposite valley, due to the larger num- 
ber of destination statesi^^ Of the / processes, the /-LA 
(phonon energy Tiwq = 46.3 me V—) and /-TO {Tiloq = 
59.1meV-2£) interactions are rapid zero-order terms in 
the scattering modeli^ Scattering rates increase rapidly 
with transition energy, and saturate above the phonon 
energy. 

The electron transfer rate from initial subband i to fi- 
nal subband / is the product of the average intersubband 
scattering rate Wi / (due to all scattering processes) and 
the initial subband population Ni. Although a simple es- 
timate of current density is proportional to the total elec- 
tron transfer rate^^^ an improved model takes account of 
the spatial separation of electrons. The current density 
is therefore 



j^!^Y,N,Y,W.f{{z)f^{z),), (16) 

^ f 



where qe is the electronic charge, Lp is the length of a 
structural period and (z) is the expectation value of the 
position operator. 

The active region gain or absorption for each transition 
was calculated as^^i 



cnoeoLp 



N,\z,f\^sgn[E,})L{u;,E,f), (17) 



where riQ is the refractive index of Si, eo is the permittiv- 
ity of free space, and c is the speed of light in vacuo. Eif 
is the energy difference between the subband minima, 
sgn(-) is the signum function, Zij is the dipole matrix el- 
ement, and L{uj,Eif) is the lineshape for the transition. 
The gain spectrum is found by summing Eqn. llTl over all 
transitions. 

Although linewidth may be obtained directly from our 
scattering rate calculations^!^' or from more sophisticated 
modelS)^ several important implementation issues, such 
as how to treat the extremely broad absorptions into 
weakly bound higher energy subbands, are beyond the 
scope of the present work. Normalized Lorentzian line- 
shapes centered about Eif have been observed in THz 
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FIG. 3: A seven-well bound-to-continuum QCL at a bias 
of 7 kVcm-i with layer widths of 3.2/3.7/0.8/6.4/1.6/5.7/ 
1.8/5.1/2.0/4.7/2.0/4.3/2.2/4.1, where boldface denotes 
40% Ge barriers and lightface denotes pure Si wells. Dopants 
are spread evenly through the structure with a concentration 
of 5 X 10^® cm~^. The conduction band potential (solid line) is 
shown, with spatially dependent probability densities super- 
imposed at each subband minimum. The upper laser subband 
is shown in bold, while other subbands are shown as dashed 
lines. 



III-V systems, with full-width at half-maximum around 
2meV at low temperatures We therefore calculated 
gain spectra using linewidths in the range 1.5-2.5 meV. 
At higher lattice temperatures, the increased scattering 
rates cause linewidth broadeningj^^*^ and consequently 
we expect our higher linewidth results to be more realis- 
tic at higher temperatures. 



V. QCL PERFORMANCE 

We have described the advantages of the (111) orienta- 
tion for Si/SiGe QCLs in general terms and have previ- 
ously predicted net gain in a novel phonon depopulation 
QCL;S In this section, we present a bound-to-continuum 
active region design with a double-metal waveguide and 
demonstrate that net gain is achievable up to 90 K. We 
also show that reasonable limitations in growth quality 
due to interdiffusion do not present a significant obstacle. 



A. Active region design 

The band structure in a seven-well bound-to- 
continuum QCL was calculated using a single band EMA 
and is shown in Fig. [5] for a 7 kVcm~^ applied electric 
field. Non-radiative depopulation of the upper laser sub- 
band was reduced by limiting the higher energy scatter- 
ing processes. Pure Si wells were used to minimize alloy 
disorder scattering, and a relatively low Ge composition 
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FIG. 4: Current density as a function of applied electric field 
for the QCL design in Fig. |3l Results are shown for linearly 
increasing lattice temperatures, in the direction of the arrow, 
between 4 and 100 K. 




2 4 6 8 

Frequency (THz) 

FIG. 5: Active region gain at 7 kVcrn"^, and 4 K lattice 
temperature as a function of THz frequency for the design in 
Fig. [3] Arrows indicate linearly increasing diffusion lengths 
from to 8 A. The resulting frequency shift of the largest 
gain peak is shown inset. 



of 40% was selected for the barriers to reduce interface 
roughness scattering. The optical transition energy was 
chosen to be significantly smaller than 46.3 meV to re- 
duce /-LA and /-TO phonon emission rates. 

As modulation doping of donors in Si/SiGe het- 
erostructures may be difficult dopants were assumed 
to be spread evenly throughout the structure. It was 
also assumed that all donors were ionized at low temper- 
atures. A donor concentration of 5 x 10^^ cm~^ (sheet 
doping density of 2.4 x 10^^ cm~^) caused negligible inter- 
nal electric fields, while still allowing rapid depopulation 
of upper miniband states by Coulombic scattering. 

Figure [5] shows gain at frequencies around 5 THz, 



FIG. 6: Peak gain near 5 THz as a function of lattice tem- 
perature for the QCL design in Fig. |3] Results are given for 
linewidths of 1.5, 2.0 and 2.5 meV. The threshold gains for 10 
and 15 fim active region thicknesses are shown as dotted and 
dashed lines respectively. 



due to transitions from the upper laser subband to up- 
per miniband states, for 4K lattice temperature and 
7kVcni^^ applied electric field. Figure S] shows that a 
current density of 2kAcm~^ corresponds to these con- 
ditions, and doubles as lattice temperature increases to 
100 K. An energy balance approach^® was used to find 
electron temperatures. Assuming an identical tempera- 
ture in each subband, we found that electron tempera- 
tures increased from 110 to 150 K as lattice temperatures 
increased from 4 to 100 K. 

Figure [S] shows the peak gain near 5 THz as a func- 
tion of lattice temperature for 1.5, 2.0 and 2.5 meV 
transition linewidths. As discussed in section IIVI 
linewidth increases with temperature, and we expect 
the low-temperature gain to be enhanced and the high- 
temperature gain to be limited accordingly. 

To provide a fair comparison with subbands in 
the (001) orientation, an equivalent (001) design was 
developed. Firstly, the conduction band offset was 
matched to that of the (111) design by setting the bar- 
rier Ge fraction to 26%. Noting that the subband ener- 
gies are much smaller than the barrier energy, an infinite 
quantum well approximation was used to transform the 
layer widths and preserve subband separation, 



,(001) 



,(111) 



(111) 

(001) ~ u-o-^-^'a; • 



(18) 



This scaling also preserves the barrier transparency, i.e. 
the coupling between adjacent wells. Finally, the applied 
electric field was increased to 13.1 kVcm^^ to account for 
the reduced period length. 

The A^ subbands in the (001) system were at suffi- 
ciently low energy compared with the A^y valleys, for 
/-phonon interactions to be negligible. With the larger 
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electric field, and with only slow (7-phonon interactions 
available to cool the electron distribution, the steady- 
state electron temperature was greatly increased to 290 K 
for a lattice temperature of 4K. Consequently, thermal 
backscattering led to a much lower peak gain coefficient of 
0.2 cm^^ at 6THz, corresponding to a current density of 
1.1 kAcm~^. Growth of such a structure is also expected 
to be challenging, as the minimum layer thickness was 
reduced to 4 A. 



B. Waveguide design 

A suitable waveguide was designed using a one- 
dimensional simulation. The propagation constant was 
obtained using a transfer matrix method"*^ and the com- 
plex permittivities were found using the bulk Drude 
model using Si/SiGe material parameters from Ref. |4^ . 
The active region was modelled as bulk Sio.gGeo.i (to 
match the virtual substrate), with a doping concentra- 
tion of 5 X 10^^ cm~^. 

Surface-plasmon configurations^ were found to be un- 
suitable, due to the low confinement factor F and large 
waveguide losses a^. We therefore chose a metal- metal 
configuration, which has proved successful in GaAs-based 
THz QCLs."*^ The active region was enclosed between a 
pair of highly doped {n = lO^^cm"^), 20 nm thick Si 
layers, followed by the metallic layers. The optical prop- 
erties of the metallic layers are given in Ref. 44. 

An initial design, using gold metallic layers and 
a 10 /im thick active region gave aw=50.7cm^^ and 
r=0.99. Assuming mirror losses of — Icm"^,**^ the 
threshold gain was determined as (7Th — (ow + am)/r = 
51.7 cm~^, which was too high to achieve lasing. 

The highest temperature operation of a GaAs-based 
THz QCL was achieved recently by using copper instead 
of gold layers.— Incorporating this into our waveguide 
reduces the threshold gain to 40.5 cm~^. As shown in 
Fig. [SI this permits lasing up to T=58K for a 2meV 
linewidth. By increasing the active region thickness 
to 15 /im, the threshold gain was reduced further to 
36.9 cm""'^, which permits lasing up to T— 90K. 

C. Growth variations 

We restricted our (111) design to layer thicknesses 
above 8 A as Si/SiGe epitaxy is not as well established 
as in III-V systems. The thinnest barrier would ideally 
be thinner than this to increase the dipole matrix ele- 
ment between the bound subband and the upper mini- 
band states. The requirement for thin barriers is less 
important in (111) heterostructures than in (001) how- 



ever, as the quantization effective mass is smaller and the 
matrix element is larger. 

Ge surface segregation has been observed in (001) 
heterostructures,'*" and presumably this will also be the 
case in (111) systems. The geometry of the thinnest lay- 
ers in a QCL is expected to change considerably as a re- 
sult. We have previously shown however, that a limited 
amount of interdiffusion is tolerable, although changes in 
transition energies are expected.— 

Figure [5] shows that gain increases slightly as a func- 
tion of interdiffusion length, Ld up to 8 A for our QCL 
design. This is due to the thinnest barrier being de- 
graded, increasing the dipole matrix element for optical 
transitions. The inset in the figure shows that the cen- 
ter of the largest gain peak correspondingly shifts from 
around 5 to 4.7 THz, as the upper laser subband energy 
decreases 1^ 



VI. CONCLUSION 

We have shown that intersubband lasing in the A val- 
leys of Si/SiGe heterostructures becomes viable in the 
(111) orientation. Although the conduction band off- 
set in the (001) orientation is large, the usable energy 
range was shown to be superior in the (111) orientation 
due to the sixfold valley degeneracy. The quantization 
effective mass was also shown to be much smaller in 
(111) heterostructures, and complications due to inter- 
valley mixing and uniaxial strain splitting are avoided. 

We have presented a bound-to-continuum design for 
a (111) Si/SiGe QCL and investigated several options 
for waveguides. A surface-plasmon waveguide was shown 
to be inadequate, while good results were achieved 
for a double-metal configuration using copper metal- 
lic layers. We have shown using a self-consistent rate- 
equation/energy balance calculation that net gain at 
5 THz is possible up to a lattice temperature of 90 K, 
with a low-temperature current density of 2kAcm^* for 
a 15 fim thick active region. 

The (111) design was found to be vastly superior to 
a (001) oriented equivalent, due to the phonon-mediated 
electron cooling and the reduced effective mass. It was 
also shown to be tolerant to, and indeed improve slightly, 
with Ge interdiffusion lengths up to 8 A. 
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